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Abstract 

The two dark sectors of the universe - dark matter and dark energy - may 
interact with each other. Background and linear density perturbation evolution 
equations are developed for a generic coupling. We then establish the general 
conditions necessary to obtain models free from early time non-adiabatic insta- 
bilities. As an application, we consider a viable universe in which the interaction 
strength is proportional to the dark energy density. The scenario does not ex- 
hibit "phantom crossing" and is free from instabilities, including early ones. A 
sizeable interaction strength is compatible with combined WMAP, HST, SN, LSS 
and H(z) data. Neutrino mass and/or cosmic curvature are allowed to be larger 
than in non-interacting models. Our analysis sheds light as well on unstable 
scenarios previously proposed. 
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1 Introduction 



Current cosmological measurements point to a flat universe whose mass-energy com- 
position includes 5% ordinary matter and 22% non-baryonic dark matter, while it 
is dominated by the so-called "dark energy" component, identified as the engine for 
accelerated expansion [1-5]. 

The most economical description of current cosmological measurements explains 
the nature of dark energy as a Cosmological Constant (CC) in Einstein's equations, 
representing an invariable vacuum energy density. The equation of state w of the dark 
energy component in the CC case is constant and w = P^j pd e = —1, where Pd e and p de 
denote dark energy pressure and density, respectively. However, when computing the 
vacuum energy density from the quantum field theory approach, the naively expected 
value exceeds the measured one by 123 orders of magnitude and it needs to be cancelled 
by extreme fine-tuning. This unhappy situation has been dubbed the CC Problem. 
Disregarding anthropic justifications, and in the absence of a fundamental symmetry 
which sets the vacuum energy to vanishingly small values, it is appropriate to look 
for alternative physical mechanisms. A related problem is the so called why now? or 
coincidence problem, i.e. why the dark matter and dark energy contributions to the 
energy budget of the universe are similar at this precise moment of the cosmic history. 

An a priori appealing avenue is to look for a dynamical explanation of the accel- 
erated expansion. Under the inspiration of the idea of inflation, it has become quite 
popular to consider cosmic scalar fields, dubbed quintessence, which would drive the 
expansion of the universe [6-11]. However, quintessence models are not better than 
the CC scenario as regards fine-tuning: no symmetry explains the tiny value of the 
potential at its minimum, which is imposed by hand. In spite of this, it seems worth to 
study the role that scalar fields may play in the evolution of the universe and explore 
the possibility of a dynamic understanding of the problem. 

Cosmic scalar fields, if present, may couple to all other fields in nature. While the 
strength of interactions between ordinary matter and the dark energy fields are severely 
constrained by observation [12], significant interactions within the dark sectors itself, 
i. e. between dark matter and dark energy, are still allowed and could affect significantly 
the universe evolution. Interacting dark matter-dark energy models have been proposed 
since the 90's (see [13] for a complete set of references). They were first explored in 
the context of coupled quintessence [9,14,15]. Some quintessence models gain extra 
motivation for being particular cases of theories of modified gravity (i.e. Brans-Dicke 
theories). Others have been proposed as a natural explanation to the coincidence 
problem [16]. Finally, it was pointed out [17] that when dark matter and dark energy 
interact, the system may mimic an effective w < — 1 naturally: current data still allows 
for such a possibility. 

In this paper we explore the simple idea that the densities of dark matter and dark 
energy do not evolve independently but coupled, although we will not refer to any 
particular cosmic field. The interaction strength between the two dark sectors will be 
generically dubbed as "dark coupling". To fix the ideas, the possible interactions can 
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be parametrized [18] by 



V M Tg m)i/ = Qu^/a, (1) 
V^d), = -Qui^/a, (2) 

with Tfc m s and T^ e ^ the energy momentum tensors for the dark matter and dark 

energy components, respectively. The dark matter four velocity iii™^ is defined in the 
synchronous gauge, in terms of the fluid proper velocity v^ dm y as uh™^ = a(—l,v dm ), 
where [i = 0..3 and i = 1..3. The coefficient Q encodes the dark coupling, that is, it 
stands generically for the interaction rate between the two dark sectors 00. 

In uncoupled models, any value is tantamount to dynamical dark energy 

and in consequence dark energy density perturbations develop, contrary to the case 
of pure vacuum energy (w = —1). The dark interaction, when present, will impinge 
on the density perturbations. Recently, using the results of Ref. [20,21], a consistent 
treatment of perturbations in models of dark coupling has been proposed. It was 
pointed out in Ref. [18] that early time instabilities may arise for constant equation 
of state, to / -1, driven by the coupling terms appearing in the non-adiabatic dark 
energy pressure perturbations. It was also claimed that such instabilities are present 
no matter how weak is the coupling. 

In this paper: 

• We formalize the evolution equations for a general coupling Q, up to first order 
in linear perturbation theory. 

• We definitely clarify the origin of the non-adiabatic instabilities, identify the 
instability regions as a function of the model parameters and propose the general 
conditions necessary to avoid them. 

• The above results will be then illustrated within a successful class of models, in 
which Q is proportional to the dark energy density. Present data will be shown to 
allow for a sizeable interaction strength and to imply weaker cosmological limits 
on neutrino masses or cosmic curvature with respect to non-interacting scenarios. 

• Finally, a comparative discussion of the existing literature will shed light on 
previously proposed models. 



6 It is to be noticed that an alternative parametrization, in which the right-hand side of Eqs. |T]) 
and ([2]) would be substituted simply by Q and — Q, would lead for the choices of Q considered in this 
work to the same results than those found below, as the difference only shows up in the equations for 
the 6dm evolution, which are not operative in the synchronous gauge, comoving with dark matter. 

7 Notice as well that more covariantly-written parametrizations are possible and even desirable, 
such as for instance to assume instead a dark coupling term oc Tjt^^v^', or oc Th^u^, with i,j 
denoting some species, and alike phenomenological ansatzs. The qualitative findings to be developed 
in this paper will also hold for such constructions, whose details will be explored elsewhere [19]. In 
order to compare with previous work in the literature we will stick in this paper to the formulation 
adopted in Eqs. (TT|) and (|2|) above. 
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2 Background and linear perturbations 



Consider a flat universe described by the Friedman-Robertson- Walker (FRW) metric. 
We will work in the synchronous gauge [22,23] comoving with dark matter, i.e. a gauge 
in which the dark matter peculiar velocity vanishes. Focusing on the dark energy and 
dark matter evolution and assuming pressureless dark matter w dm = Pdm/ Pdm = 0, it 
follows from Eqs. ([1]) and (j2J) that the evolution of the background energy densities is 
given by: 

Pdm + %Hpdm — Q j (3) 
pde + 3n Pde {l + w) = -Q, (4) 

where the dot indicates derivative with respect to conformal time dr = dt/a and 
H = a/a. The effective background equations of state for the two fluids are thus given 
by 

W dm = -3&: ' ( 5 ) 

"S'^ + sfc. (6) 
where w would be the dark-energy equation of state for a vanishing interaction. Equa- 
tions ([5]) and ([6]) suggest immediately how the interaction between the two fluids can 
contribute to the effective value of the dark energy equation of state. Furthermore, they 
show that positive (negative) values of Q contribute as an effective negative (positive) 
pressure in the dark matter background equation. This leads to less (more) dark mat- 
ter in the past than in the uncoupled case and, as a consequence, the matter radiation 
equality will happen later (earlier) on. 

Also notice that a universe in accelerated expansion today requires w < —1/3 even 
in the presence of a dark coupling. Indeed, the deceleration parameter satisfies: 

q = -^ = l(l + 3wn de ) (7) 

either with or without dark coupling. In this equation we neglected the curvature 
contribution. 

The coupling between the two dark sectors will also affect the evolution of the 
dark matter and dark energy density perturbations, Sp dm and Sp de , respectively. In 
the synchronous comoving gauge, metric scalar perturbations are described by the two 
usual fields [24] h(x,r) and t](x,t). Defining 5 = Sp/p for the fluid density perturba- 
tions, 9 = diV 1 for the divergence of the fluid proper velocity v 1 and using Eq. ([!]) for 
presureless dark matter, it results at first order in perturbation theory: 

1 . 

<W = --h + S[Q/p dm ] . (8) 
For dark energy, using Eq. ([2]) it follows that: 

S de = -(l + w)(e de + ^h)-3n(^-wjS de -5[Q/p de \ , (9) 

e de = -n(i - 3w)e de + - k " 6 f de 5 de + ^-e de . (10) 

1 + W dp de p de 
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The initial conditions for the dark matter and dark energy density perturbations will 
be taken in what follows as in Ref [24]; in particular we assume 5 de = initially. For 
numerical computations and illustrative plots, we have used the publicly available CAMB 
code [25], which includes the full evolution for all species (photons, neutrinos, baryons, 
dark matter and dark energy), modifying it to take into account the dark coupling. 



2.1 The doom factor 

Equations (jSj), (Q and (flOj) are written in the dark matter rest frame. However, a 

SPde 
&Pde 



priori, the dark energy sound speed, c 2 sde = j-^, is only well known in the rest frame 



of dark energy [26]. It can be shown [26] that 

sp de = c 2 sde Sp, 



Me 



_ 2 \ • ®de 

sde C a del Pde ^ 2 ' 



'ID 



where c 2 de is the propagation speed of pressure fluctuations in the rest frame of dark 
energy and c 2 ade = P de / p de is the so called "adiabatic sound speed", which for constant 
w satisfies c 2 ade = w. In the presence of a dark coupling, from Eqs. (0J and (TTTj) it 
follows that 



SP, 



dc 



$Pde 



+ 3(c 2 sde-C 2 ade )(l + W e // 



no, 



de 



k 2 5, 



de 



where we define 



C 2 sde + 3{c 2 sde -C 2 ade ){l+ W ) (1 + d; 



Q 



ne, 



de 



k 2 5, 



de 



(12) 



(13) 



3Hp de {l + w) 

where d stands for doom: we dub it so as it is precisely this extra factor, proportional 
to the dark coupling Q, which may induce non-adiabatic instabilities in the evolution 
of dark energy perturbations. Its sign will be determinant, as we are going to show. 



Rewriting Eqs. ([8]), (Q and ffTUl) in terms of c 2 de and d, we have: 



$dm 
Sde 



-~h + 3H(l + w)5 

2 \_pdn 

-(l + w)(6 de + h)-3H(l + w)5[d} 



(14) 



-3W {c 2 de - w) 



6 de + 3H(1 + w) (1 + d) 



7 de 



-H (l - 3c 2 sde - 3d(c 2 sde + 1)) 6 de + — — c 2 de 5 de 

' 1 + w 



(15) 
(16) 



where S[d] includes 5Q and S de contributions. 

Below, in addition to w < — 1/3, the speed of sound c 2 de will be assumed positive, 
with c 2 de = 1 in numerical computations. 
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2.2 Early time (in)stabilities 



In general, the evolution of the dark energy and the dark matter perturbations are 
directly coupled. Indeed, it is known that, even in the uncoupled case, once w/-l 
the dark matter and dark energy perturbation evolution depend on each other [27] (see 
also Sec. 12. 3j) . 

When the dark fluid linear perturbation equations are combined into second order 
differential equations, they take the generic form [5 

b" dm = A m % + B m & + F{p u 5 h 5'-, i ± dm) , (17) 
cr a 

S'L = + B e S -^ + g(p i ,5 i ,5' i ;i^de), (18) 

a z a 

where ' = d/da and the function T {Q) stores the dependence in all variables but 8dm 
or $dm ( S de or 8' de ). 

The evolution of a perturbation will depend on the relative weight of the three 
terms in the corresponding equation and on their signs: 

1. For positive A, the A and B terms taken by themselves would induce a rapid 
growth of the perturbation, which may be damped or antidamped (reinforced) 
depending on whether B is negative or positive, respectively |^. In particular, 
for A and B both positive, the solution may enter in an exponentially growing, 
unstable, regime. 

2. For negative A, in contrast, the A and B terms taken alone describe a harmonic 
oscillator, with oscillations damped (antidamped) if B is negative (positive). In 
the A, B < regime, the third term may plays in fact the leading role. 

It is worth reviewing the uncoupled scenario in detail before proceeding further: in it, 
dark matter perturbations behave as in case [1] above (with A > and B < 0), while 
dark energy ones provide an example of behavior as in case [2j 



2.3 Uncoupled case 

Consider first the growth of dark matter and dark energy density perturbations, in the 
absence of dark coupling (Q = 0), at large scales (H/k 2 <C 1) and early times, when 
Qde^de can be neglected as Vtde -C fli for i = all other species. In this regime, 

x" ^ dm ^ 5 d . . 

°dm = n S W-r - 9— + ^' (I 9 ) 

La CL _ ( ' 

s'L = -l(cL e -w) 6 -^-( 5 --3w) 6 f + g, (20) 

8 We surrender here to the extend habit of expressing first order temporal differential equations in 
terms of conformal time, ' = d/dr, while in second order ones ' — d/da is used. 

9 Obviously, for \B\ >> \A\ a negative B would prevent the onset of growth for any sign of A. 
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Figure 1: Left panel: The blue (red) curve depicts the evolution of the Sd m perturbation 
vs the scale factor for k = 0.001 h/Mpc and w = —1.1 (w = —0.9). Right panel: same 
as in the left panel but for the 5d e evolution. 



with the functions T and Q are given by: 



T 



[1 + W 



£ 1+3 



8 Pi 



Si 



3E(l+3g)4-2(3 S; 



de 



ft 

°dm 



_ i=£de 



(21) 



Equation (fl9j) for dark matter density perturbations has A > and B < 0: the 
latter term damps the growth of dark matter perturbations propitiated by the A term, 
with the overall well known polynomial rising, see Fig. [1] (left panel) . 

Equation (1201 for dark energy perturbations has instead both A and B coefficients 
negative, with the first two terms describing then a damped harmonic oscillator and 
in this case the contribution of radiation, dark matter and matter encoded in Q drives 
the evolution, see Fig. [1] (right panel). Equations ff20|) and fl2TT) also illustrate that, for 
constant w ^ — 1, dark energy perturbations do develop even in absence of coupling, 
seeded by the Q term, in contrast to the w = — 1 (pure vacuum energy) case, in which 
no perturbation can develop in the dark energy background. 
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2.4 Strongly coupled case 



Consider now the opposite case in which the dark-coupling terms dominate over the 
usual one. For c 2 sde > considered all through, the strong coupling regime can be 
characterized by 



Q 



> 1, (22) 



which guarantees that the interaction among the two dark sectors drives the non- 
adiabatic contribution to the dark energy pressure wave, see Eq. ( TT2i) . and becomes the 
leading term in Eq. ( 1T51) for 5d e and Eq. ( |T6l) for Ode- At large scales, those equations 
reduce to: 

5 de ~ -3H(c 2 sde -w) (s de + 3H(l + w)d^j -(1 + w) ^3W<J[d] + |j ,(23) 

de ~ 3Hd(cL e + l)^ e + ^^^ e . (24) 

The complete second order differential equation describing the growth of dark energy 
perturbations can be found in Appendix |Aj for values of w near —1 and in the strong 
coupling regime, the Sd e and S de contribution to the second order differential equation 
reads 

<4 , o M3«fc_-«0 , 3(1 + w) 

2 

's de 



b" de ~ 3d(cL e + l)( ^ + 3^ v y ^ +^^[d] )+... (25) 



The sign of the coefficient B e of S' de in this expression is crucial for the analysis of 
instabilities, as previously argued, see Eq. ([IB]) . Notice that the size and sign of the 
B e coefficient also determines the growth rate of e in Eq. ([241 . Assuming c^ de > 0, it 
reduces to the sign of the doom factor d defined in Eq. ( TT31 . 

Indeed, as previously argued, a positive d acts as an antidamping source in the 
growth Eq. (1251) . Whenever d > 1, it will trigger an exponential runaway growth of the 
dark energy perturbations when simultaneously the overall sign of the A e coefficient 
of d~ de , resulting from the last two terms in Eq. fT25|) . is also positive. Large scale 
instabilities arise then and the universe appears to be nonviable. 



3 A simple viable model: Q oc pd e 

Let us consider now a specific simple coupled model, with the interaction rate Q = 
^TCpde- Equations (CQ) and (j2D become consequently: 

V^ m) „ = iHpdeU^/a, (26) 
VmZJU, = -in Pde u^/a, (27) 
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where £ is a dimensionless coupling which will be taken as constant. It parametrizes 
in this model the dark coupling strength. Somewhat similar coupled models have been 
explored in the literature [15,18,28-31] In particular the authors of Ref. [18], assuming 
Q Pdm or Q oc pdm + Pde, concluded that their models, with constant w, do not 
provide viable scenarios due to the presence of early time instabilities. 

We will show here, though, that the model described by Eqs. (126]) and (1271) can 
satisfy all current observational constraints, without suffering from the instabilities 
pointed out in [18], even for constant w. In this section we will analyze this model in 
detail, postponing, and extending, the discussion of the scenarios presented in Ref. [18] 
to Sec. H 



3.1 Background 

The dark matter and dark energy background densities evolve in this model as two 
fluids, coupled with a strength linearly dependent on the dark energy density present 
at any given time of the cosmic evolution, 

Pdm + 3Hp dm = iT-i-Pde (28) 
f) de + 3Hpde(l + Wde) = -£,Ti-Pde, (29) 

which then correspond to two fluids with effective equations of state given by 

= ( 30 ) 

w e JJ = w + l. (31) 

Note that for constant w, w e JJ is also constant while w e J^ is redshift dependent. The 
solutions to equations ([3]) and (j4]) are then: 



Pdm = PZ^ + P^-^jil-a-^)^ 6 , (32) 

iW de 

P d e = Pf e a-^JJ\ (33) 



The dark energy density is thus always positive, all along the cosmic evolution and 
since its initial moment. To ensure that the same happens with the dark matter density, 
all values of w < are acceptable for £ < 0, while for positive £ it is required that 

The resulting evolution of the relative energy densities is shown in Fig. [21 where 
the present matter and dark energy densities have been imposed to be 0.27 and 0.73 
(a procedure to be repeated all through the paper). We see that negative (positive) 
couplings lead to more (less) dark matter in the past than in the uncoupled case, as 
expected. As a consequence, the matter radiation equality happens earlier (later) on. 
Fig. (j2j) illustrates as well that positive (negative) values of Q soften (worsen) the 
coincidence problem, although the difference turns out to be quantitatively minor for 
the phenomenologically allowed values of £. It will be shown in Sect. 13.51 that models 
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a 



Figure 2: Scenario with Q oc pd e - Relative energy densities of dark matter plus baryons 
Qdm+b (blue) , radiation Q ra d (black) and dark energy Q de (red), as a function of the 
scale factor a, for w=-0.9. Three values of the coupling are illustrated: £ = (solid 
curve), 0.25 (long dashed curve) and —0.25 (short dashed curve). 



with £ < and 1 + w > will give the best agreement with large scale structures 
observations. 

Finally, Eqs. fl3T|) and fl33l) show that the so-called phantom regime, in which the 
dark energy density would diverge in the future, can only happen in this model for 
1 + w e JJ < 0, that is, for w < — 1 — £/3, instead of w < — 1 in the uncoupled case. 

3.2 Dark coupling vs. dynamical dark energy 

From the previous analysis one can wonder what would be the dark energy equation of 
state reconstructed from observational data, if analyzed assuming no coupling. Indeed 
a dynamical, redshift-dependent, equation of state w(z) (where z denotes redshift), can 
be mimicked by the combination of constant w plus the dark coupling £ . 

A wealth of data are sensitive to the Hubble parameter H(z) or functions of its inte- 
gral. When such data are analyzed assuming no dark coupling, the following expression 
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for the Friedmann equation is to be considered q 



HHz 



Rh{z) = = n£> (1 + z )» + exp 



3 / dz 



,1 + w{z') 
1 + z' 



(34) 



where flf^ = pf / Pc\ being the critical energy density today. Relation (134"1) can 
be inverted to obtain w(z) as a function of Rh(z): 

_ 1 (1 + z)dR H /dz - 3R H 

() "3 ^-n£(i + ,)3 • (35) 

An analogous relation between the redshift dependent equation of state and the lumi- 
nosity distance and its derivative was also obtained previously [32]. 

In contrast, in the presence of the dark coupling, from Eqs. (j32l and fl33l) it results 
for the simple model analyzed in this section 



R H (z) = (1 + z 



,3 



n£ + ^Aj (i - (i + *) 3w ") + + *)■■*' 



.(36) 



Comparing Eqs. ( |34|) and ( 1361) . it follows the relation between the hypothetically re- 
constructed expression of a dynamical equation of state w(z) and the constant w,£ 
parameters of the darkly coupled universe: 

™( z ) = e — : — : o , (37) 

1- -4tt (l-(l + « - 3tt * 

an expression which at small redshifts tends to: 

w(z) ~ w(l + £z) , (38) 

while w(^) ^0 at very large redshifts. 

A striking implication of Eq. fl37|) is that the reconstructed equation of state would 
show a peculiar divergent behavior for those models in which the denominator vanishes. 
In fact, for modified gravity models and coupled scalar-tensor models a similar behavior 
was previously pointed out [33,34]: their reconstructed w(z) could appear singular at 
some point in the past and cross the "phantom boundary" w = —1. 

In the phenomenologically viable model discussed here no such divergent behavior 
can arise in the w(z) reconstruction, for negative values of the dark coupling £ for 
any w < 0. Whereas, for positive dark coupling values, a "phantom-like" behavior 
appears, with the divergence occurring at redshift z ~ (9(1 — 10), depending on the 
specific values of the parameters. 

All these different behaviors are illustrated in Fig. [31 for w = —0.9 and £ = 
±0.8, ±0.2. The model with w = —0.9 and £ = 0.8 is not viable, due to negative 
dark matter energy density, and it is shown only for illustration purposes. Notice that 



10 We obviate here as well, for the sake of the argument, the contribution to Q^ot °f an y species 
other than dark matter and dark energy: baryons, neutrinos and curvature. 
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Figure 3: Scenario with Q oc pd e - Reconstructed w(z) as function of z, for w = —0.9. 
The black (solid) and magenta (short dashed) curves depict the w(z) behaviour for 
£ = 0.2 and £ = 0.8. The blue (long-short dashed) and the red (long dashed) curves 
denote the w(z) behaviour for t; = —0.2 and £ = —0.8. Notice that for positive values 
of £ we recover the divergent phantom- crossing behavior appearing in scalar-tensor 
theories. 



an apparently diverging behavior of the equation of state could thus be a clean indi- 
cator of the sign of a putative dark coupling. Furthermore, the degeneracy between 
a dynamical equation of state and a constant w plus dark coupling, discussed in this 
subsection, can be disentangled with other z— sensitive data, such as precision BAO 
or SuperNovae measurements and the growth of perturbations seeding large structure 
formation, as it will be seen in Sec. 13.51 

3.3 Linear perturbation theory 

The propagation of dark energy pressure waves is driven by Eq. (fl2j) as explained 
above, with the doom factor Eq. ({TBI given in this case by: 



The evolution of the perturbations, in the synchronous comoving gauge, is described 
by Eqs. (fT4"|) . ffTBT) and ffTBT) with, in the present scenario, 



d 




(39) 



S[d)=0, 



(40) 




(41) 
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The dark matter density perturbation changes in size with respect to the uncoupled 
case. It is easy to confirm numerically that the major effect of the coupling on the dark 
matter distribution results from the background, though, for dark coupling values small 
enough so as to allow viable models. For negative couplings, there is more dark matter 
in the past. As a consequence, dark matter perturbations cluster more and the density 
perturbation is larger. 

3.4 Early time (in)stability 

The strong coupling regime is defined by |d| > 1, as explained in Sec. 12.41 In this 
regime and for large scales, 7i/k ^> 1, the dark energy perturbations are described by 
Eqs. (1231) and (!24j) . where in this case Eqs. (I39I) . (140]) and (l4"T!) above apply. It leads, 
for constant w ^ — 1, to a 8d e and S de contribution to the early time growth at large 
scales, described by 

6 de * 3(1 + c sde ) d ^ - + 3-,^-^ J + ... (42) 

which shows that the sign of d defines the (un) stable regimes ["1: 

1. For d < 0, that is, for £ < and 1 + u; > (or £ >0 and 1 + w < 0), no 
instabilities are expected; 

2. When £ and 1 + w have the same sign, instabilities can develop at early times 
whenever d > 1 . 

Sweeping over all possible values of £ and w, we have indeed confirmed numerically 
these predictions. We summarize our instability criteria for the present model in Tab. [TJ 
An example of the onset of unstable behaviour is shown in Fig. HI right panel, in which 
the particular choice of w and £ makes the model unstable at early times, see Tab. 1. 
The left panel of Fig. @] depicts instead the dark energy perturbation, 5d e , versus the 
scale factor a in a model free of instabilities and pathological behaviors in the linear 
perturbation evolution, notwithstanding the strong coupling regime. 

In the next subsection, we will confront in detail the model with data, sweeping 
over values of the dark coupling strength £ < 0, although models with £ > are also 
viable as long as £ > —w, as previously shown. 

3.5 Cosmological constraints 

We explore here the current constraints on the dark energy-dark matter coupling £, 
considered constant, allowing for a non zero spatial curvature. The framework used is 
a cosmological model described by ten free parameters, 

6 = {U b , LO dm , 9c MB: fifc, fu, w, €, n s, A s } , (43) 
11 Recall that c 2 sde > is assumed all through. 
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Model: Q oc pd e 


1 + w 




Pdm 


Pde 


d 


Early time 
instability? 




+ 


+ 


T 


+ 


+ 


Yes 




+ 




+ 


+ 




No 








+ 


+ 


+ 


Yes 






+ 


T 


+ 




No 



Table 1: Scenario with Q oc pd e - Stability criteria driven by the sign of d in Eq. ( T?Pj) . 
for |d| > 1. The =p szgns indicates that pd m is negative in the past for large positive 
couplings only. 



Parameter 


Prior 




0.005-0.1 


^dm 


0.01-0.99 


6 CM B 


0.5-10 


T 


0.01-0.8 


Q k 


-0.1-0.1 


fu 


0-0.3 


W 


-1-0 




-2-0 


n s 


0.5-1.5 


ln(10 10 A s ) 


2.7-4.0 



Table 2: Priors for the cosmological fit parameters considered in this work. All priors 
are uniform in the given intervals. 



where CO}, = Qf,h 2 and uj,im = Qdmh 2 are the physical baryon and dark matter densities 
respectively, #cm£0 is proportional to the ratio of the sound horizon to the angular 
diameter distance, r is the reionisation optical depth, Q k is the spatial curvature, 
/„ = Q u /Qdm refers to the neutrino fraction, n s is the scalar spectral index and A s the 
scalar amplitude. The priors adopted on those parameters are given in Tab. [2j 

We use the publicly available package cosmomc [35], modifying it in order to in- 
clude the coupling among the dark matter and dark energy components, for the model 
considered in this section. 

A conservative compendium of cosmological datasets is considered. First, what we 
call run includes the WMAP 5-year data [1,2], a prior on the Hubble parameter 
of 72 ± 8 km/s/Mpc from the Hubble key project (HST) [36], the constraints coming 
from the latest compilation of supernovae [3] and the H(z) dataEf] at < z < 1.8 from 
galaxy ages obtained by [37]. We present a second data analysis, that we call run I, in 

12 The 9c mb parameter can be replaced by the Hq parameter. However, using 8cm b is better due 
to its smaller correlation with the remaining parameters. 

13 We thank R. Jimenez and L. Verde for suggesting the H(z) data addition to the analysis. 



13 




a, scale factor a, scale factor 



Figure 4: Scenario with Q oc pd e in the strong coupling regime. Left (right) panel: 
Evolution of the 5d e perturbation vs the scale factor for scales k = 0.001 h/Mpc, for 
£ = —0.8 and w = —0.9 (w = —1.1 ). Notice the early time instability present when the 
doom factor d in Eq. / fffPj) is sizeable and positive, as predicted from the study of the 
instabilities. 

which we add to run the data on the matter power spectrum (large scale structure 
data or LSS data) from the spectroscopic survey of Luminous Red Galaxies (LRGs) 
from the Sloan Digital Sky Survey (SDSS) survey [4]. 
In summary, 

• nm0=WMAP(5yr)+HST+SN +H(z), 

• runl=run0 +LSS. 

Figure (left panel) illustrates the 1 and 2a marginalized contours in the £,-Qdmh 2 
plane. The results from the two runs described above are shown. We restrict ourselves 
here to negative couplings and w > — 1, which guarantees that instability problems in 
the dark energy perturbation equations are avoided for all values of £. Notice that a 
huge degeneracy is present, being £ and fldmh 2 positively correlated. The shape of the 
contours can be easily understood. In a universe with a negative dark coupling £, the 
matter content in the past is higher than in the standard ACDM scenario due to an 
extra contribution proportional to the dark energy component, see Eq. ( I32p . Therefore, 
the amount of intrinsic dark matter needed - that is, not including the contribution of 
dark energy through the coupling term - should decrease as the dark coupling becomes 
more and more negative and can be as small as 0.02, as indicated by run results. The 
addition of LSS data to the analysis reduces considerably the allowed parameter space. 
This is due to the enormous growth of clustering for values of the coupling £ < —0.5. 
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Figure 5: Scenario with Q oc p de - Left (right) panel: la and 2a marginalized contours 
in the t;-Q dm h 2 (£-&k) plane. The largest, green contours show the current constraints 
from WMAP (5 year data), HST, SN and H(z) data. The smallest, red contours show 
the current constraints from WMAP (5 year data), HST, SN, H(z) and LSS data. 



For this range, the amplitude of the fluctuations increases, reaching values of as > 2 
and therefore providing a bad fit to LSS data0. 

The right panel of Fig. [5] shows a positive correlation between the coupling £ and 
the spatial curvature Qk- High precision CMB data indicates that currently the spatial 
curvature is a subdominant contribution to the energy budget of the universe^, which 
implies + f2$^ + £1^ ~ 1. A non zero spatial curvature component implies instead 
^de + ^dm + + ^fc = 1- A negative coupling £ will increase the dark matter 
contribution and therefore a small negative curvature (closed universe) is needed to 
compensate the effect and describe well CMB data. The degeneracy between £ and flk 
gets alleviated if one adds LSS data to the analysis. 

Figure E] (left panel) depicts the constraints on the £-w plane. We restrict ourselves 
here to w > — 1 and £ < 0, a parameter region which ensures a negative doom factor, 
see Eq. (139|) . and thus spans an instability-free region of scenarios to explore (see 
Sec. 13.41) . Current data is unable to set strong constraints on the equation of state 
parameter w. 

The right panel of Fig. O shows next the correlation among the fraction of matter 
energy-density in the form of massive neutrinos f v and the dark coupling £. The relation 
between the neutrino fraction used here f v and the neutrino mass for iVj, degenerate 
neutrinos reads 

tt u h 2 Y, m v 1 N v m v 1 , s 

U ~ n dm h 2 ~ 93.2eV ' n dm h 2 ~ 93.2eV ' tt dm h 2 ' 1 ' 

14 The authors of Ref. [2] have reported a best fit value as = 0.812 ± 0.026. 

15 The authors of Ref. [38] have found f^ 0) = -0.002±°;°^ (95% CL limits) assuming a dynamical 
dark energy component. 
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Figure 6: Scenario with Q oc pd e - Left (right) panel: la and 2a marginalized contours 
in the (£,-fv) plane. The largest, green contours show the current constraints from 
WMAP (5 year data), HST, SN and H(z) data. The smallest, red contours show the 
current constraints from WMAP (5 year data), HST, SN, H(z) and LSS data. 



Neutrinos can indeed play a relevant role in large scale structure formation and leave 
key signatures in several cosmological data sets, see Ref. [39] and references therein. 
More specifically, the amount of primordial relativistic neutrinos changes the epoch of 
the matter-radiation equality, leaving an imprint on both CMB anisotropics (through 
the so-called Integrated Sachs-Wolfe effect) and on structure formation, while non- 
relativistic neutrinos in the recent Universe suppress the growth of matter density 
fluctuations and galaxy clustering. This can be observed in Fig. where the dotted 
curve depicts the matter power spectrum for three degenerate massive neutrinos (N u = 
3) with m v = 0.4 eV. Notice that the matter power spectrum is reduced with respect 
to the m u = case, especially after the matter-radiation equality era (imprinted in 
the power spectrum as a turnover). 

There is a strong and very well known degeneracy in the ^2 m v — w plane, as first 
noticed in Ref. [40]. Cosmological neutrino mass bounds become weaker if the dark 
energy equation of state is taken as a free parameter. If w is allowed to vary, Qdm 
can take very high values, as required when m v is increased in order to have the same 
matter power spectrum. More recently, the authors of Ref. [41] have pointed out that a 
higher neutrino mass is possible if dark matter and dark energy are coupled. Figure [7] 
shows the matter power spectrum in several scenarios. The long dashed line refers to 
a universe with three massless neutrinos and a coupling in the dark sector £ = —0.24. 
Notice that the matter power spectrum is enhanced with respect to the £ = scenario, 
due to the higher matter energy density in models with a non-negligible coupling. 
However, the power enhancement effect induced by the presence of a coupling can be 
compensated by adding massive neutrinos in the game. Those neutrinos will reduce the 
power spectrum, see the short dashed curve in Fig. which is indistinguishable from 
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0.0001 0.001 0.01 0.1 1 10 

k[hMpc"'] 

Figure 7: Scenario with Q oc pd e - The black solid line depicts the matter power spectrum 
for the ACDM model with three massless neutrinos. The blue dotted line shows the 
mater power spectrum for a ACDM universe with three massive degenerate neutrinos 
with Yl m v = 1-2 e V- The green long-dashed curve refers to a universe with three 
massless neutrinos with a dark energy-dark matter coupling £ = —0.24. The red short- 
dashed line illustrates a universe with J2m u = 1.2 eV and a coupling £ = —0.24. 



the matter power spectrum in a ACDM universe. This m v — £ degeneracy is shown 
in Fig. [6] (right panel): a neutrino mass of ^2m u ~ 1.5 eV is allowed for couplings 
£ > -0.6 at the 2a level. 

4 Comparison with the literature: Q oc p^ m 

Early time non-adiabatic instabilities were pointed out in Ref. [18], in which several 
coupled models were considered. In particular they studied an interaction rate of 
the form: Q = ^Hpdmi with constant w and £, concluding that it was unstable for 
constant w even for very small values of the coupling. We will readdress that model 
here, to clarify if and when it is subject to non-adiabatic instabilities and explore some 
supplementary aspects. 
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Figure 8: Scenario with Q oc p dm . Relative energy densities of dark matter + baryon 
Qdm+b (blue), radiation VL rad (black) and dark energy Q de (red), as a function of the 
scale factor a, for w = —0.9. Three values of the coupling are illustrated: £ = (solid 
curve), 0.25 (long dashed curve) and —0.25 (short dashed curve). 



In this scenario, Eqs. ([I]) and (j2J) lead to the following background equations: 

p dm + 3Hp dm = £Hp dm , (45) 
p de + 3Hp de (l + w) = -iHp dm - (46) 

These two fluids exhibit thus effective equations of state given by 

<% = -l (47) 

< f e f = ™ + i^- (48) 

Notice that now the effective dark energy equation of state is redshift- dependent, and 
consequently, for constant w and £, 

Pdm = p2y 3+? ; (49) 

Pd e = Pd Q) e a^ 1+w) + (a- 3 ^ - a-^) . (50) 

The authors of Ref. [18] chose to restrict their analysis to constant negative £ and 
1 + w > 0. For such parameter values, in their model the dark energy density p de 
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is negative for any w < 00 This is illustrated in Fig. (jSj), where the short dashed 
curve corresponds to £ = —0.25 values and fide is seen to be negative since the universe 
birth until late times, while fldm+b is always positive. This is in stark contrast with 
the behaviour in the scenario with Q cx pd e discussed in Sect. El for which all energy 
densities remain positively defined all through the universe history. Figure flSJ) exem- 
plifies numerically as well how positive (negative) values of Q ameliorate (worsen) the 
coincidence problem, although again the effect is quantitatively unimportant for the 
phenomenologically allowed values of £. 

Moreover, the authors of Ref [18] found that the instabilities appear no matter how 
weak the coupling is. Before turning to the analysis of the doom factor and instabilities 
for this model, let us explore some further background-dependent properties. 



4.1 Reconstructing w{z) 

For this scenario, the reconstructed equation of state w(z), if obtained from data mainly 
sensitive to the fluid background and analyzed assuming no dark coupling, will diverge 
for positive couplings and redshifts around z ~ 0(1 — 10). It will show a phantom 
crossing behavior typical of scalar-tensor dark energy models, see Fig. [91 Indeed using 
Eq. (I3"5"j) . we get a somewhat more complicated expression than for the Q cx pd e scenario: 

~ n= £-(l + z) 3w+ g [(£ + 3w)r + £] 

W{Z) W -3w + (l + zy(3w + 0-(l + z) 3w+ mZ + 3w)r + £] 1 ' 

with r = ^de /^dm- ^ sm all redshifts it reduces to: 

w(z) ~w(l + -z] , (52) 



r 

while at large redshifts the reconstructed w(z) in Eq. (i5Tj) obeys 

-£ l 

w(z)2t-±- 7 ^ for £>0, (53) 



3 1 



| for £<0, (54) 



assuming |£| < 1. 



4.2 Linear perturbation theory 

The propagation of dark energy pressure waves is again driven by Eq. ( TTSl) . with the 
doom factor Eq. ( fl~3j) reading now: 

4 Pdm (55) 



3(1 + W) Pde 



16 And thus not only within the —2/3 > w > — 1 region explored in Ref. [18]. Notice that in our 
convention £ = — a, with a defined in that work. 
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Figure 9: Scenario with Q oc pdm- Reconstructed w(z) as function of z. We have 
considered w = —0.9. The black (solid) and magenta (short dashed) curves depict the 
w(z) behaviour for C, = 0.2 and £ = 0.8. The blue (long-short dashed) and the red (long- 
dashed) curves denote the w(z) behaviour for £ = —0.2 and £ = —0.8. Notice that for 
positive values of £ we recover the divergent phantom- crossing behavior appearing in 
scalar-tensor theories. 



This expression should be compared to the analogous one when Q oc pd e , see Eq.( 
they differ in the Pdm/Pde factor in the dark coupling term, which can be large at early 
times when pd e may become negligible. At early times, the doom factor will always 
thus be large even for tiny values of £, as the ratio pdm/ pde is very large then. This 
ratio is the major difference between the scenarios with Q oc pd m (or Q oc pd m + Pde 
and alike) and those in which Q oc pde- As the true expansion factor is d instead of £, 
this fact also explains why in the present scenario the phenomenological analysis only 
allows for £ values much smaller in magnitude than those permitted for Q oc pde- 

The evolution of the perturbations, in the synchronous comoving gauge, is described 
in this scenario by Eqs. (TI4")) . (TT5l) and (|T6l) with 



Mdl 



Pdr, 



3(1 + IW) pde 



0. 



<6, 



dm 



s, 



del 



(56) 
(57) 



4.3 Early time (in)stability 

The strong coupling regime is defined by |d| > 1, as explained in Sect. 12.41 In this 
regime and for large scales, Tl/k ^> 1, the dark energy perturbations are described by 
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Model: Q oc p dm 


1 + w 




Pdm 


Pde 


d 


Early time 
instability? 




+ 




+ 


T 


+ 


Yes 




+ 


+ 


+ 


+ 


+ 


Yes 






+ 


+ 


+ 




No 








+ 






No 



Table 3: Scenario with Q oc pd m - Stability criteria driven by the sign of d in Eq. (f55j). 
whenever \d\ > 1. T/ie =p szgns indicates that pd e is negative in the past for negative 
couplings independently of their value. 



Eqs. ( 1231) and (T2"4"l) . for which in this case Eqs. (}55l) . (|56|) and (|57|) above apply. It leads, 
for constant w ^ — 1, to a 5d e and 5' de contribution to the early time growth at large 
scales, described by 



S'L ^ 3(cL P + l)d 



3 ~sde _ ( 1+w)d j + 

0„ + 1 /a 2 a J 



+ ... (56 



where dots account for all other terms, which will be subdominant whenever the two 
terms above induce by themselves an exponential growth. 

Once again, for d > the 5' de term signals an antidamping -growing- regime, which 
may induce instabilities, whenever d > 1, when combined with either a positive or 
a negligible negative 5d e coefficient. Interestingly, Eq. (|5"51) shows that the latter is 
determined in this case by a competition between the two factors within brackets. We 
have verified numerically that, for values of |£| < 1 and w around -1, both terms are 
of the same order and strongly cancel. For positive d, all cases of Sec. 14.31 with A e > 
or < and B e > develop then a strong instability at early time, see upper panels in 
Fig. [TDJ 

Turning back to the characteristics of the doom factor in this model, Eq. fl55|) . it 
is to be noted that the ratio pd m / Pde not only enhances in the past the magnitude of 
d and thus the onset of the strong coupling regime even for very small £ values, but it 
also influences its sign. Indeed, recall that pd e is negative in the past for negative £, as 
can be seen from Eq. (|5"U|) . with an obvious impact on the sign of d. 

The identification of the d factor allows thus to predict the range of parameters 
for which this model will be stable, even with constant w. We have summarized the 
analysis of the early time (in)stability criteria for a model with Q oc pd m in Tab. [31 
It is illustrative and amusing to further analyze these patterns with some examples. 
Figure [TU] shows two cases of early time instabilities in the upper panels, and two 
examples of early time stability in the lower panels: 

• The two upper graphics both assume 1 + w > 0, while the dark coupling strength 
£ has opposite sign in the right and left panels. Nevertheless, both scenarios are 
unstable from their birth, as the sign of pd e is also opposite as can be seen from 
Eq. (1501) and depicted as a discontinuous line in the figures, resulting in a doom 
factor d - Eq. (1551) - which is positive for both scenarios. 
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Figure 10: Q oc pd m - Upper panels: Evolution for the 5d e perturbation vs the scale 
factor at scales k = 10 -5 h/Mpc for w = —0.9 and £ = ±0.1. Early time instability is 
present in both cases, as predicted from our analytical study. The black dashed curve 
shows the positive (negative) value of the dark energy density flde(a) (right scale). The 
lower right (left) panel shows the Sd e perturbation vs the scale factor at scales k = 10~ 5 
(k = 0.001 ) h/Mpc forw = —1.1. For both positive and negative couplings the model is 
free of early-time instabilities. However, for negative couplings, there exists a late-time 
instability (see right panel), caused by a change in the sign of the dark energy density 
fide (a), shown by the black dashed curve (right scale). 
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• The two lower panels analyze the analogous cases although with 1+w < instead. 
Again, at early times both models have opposite sign as well for the dark energy 
density pd e , combining into a safe negative value for d. 

• A late time instability is seen to appear, though, for the model in the lower right 
panel. Indeed, this is well understood as for this model pd e is seen to change sign 
during cosmic evolution, and in consequence precisely at that temporal point the 
model finds its doomsday, as d becomes positive. 

This latter example pinpoints that our (in) stability criteria on the growth of dark 
energy perturbations also applies to late evolution, if handled with care. 

Finally, as pointed out in Ref. [18], a possible way to avoid the early time instability 
of their models is to relax the assumption of a constant equation of state of the dark 
energy component and allow instead for a dynamical behavior, i.e, for a time-varying 
w(a), even when the dark coupling is turned off. For instance, the CPL [42,43] pa- 
rameterization w(a) = Wq + w a (l — a) can alleviate the instability problem because 
of its smaller ratio Pdm/Pde in the past. From our perspective, this may allow those 
models to avoid in the past the strong coupling regime as we defined it, trading it for a 
softly coupled model, that is, one in which the evolution is dominated by the coupling- 
independent terms and thus alike to the stable uncoupled regime discussed in Sec. 12.31 
A detailed perturbation analysis with different dark energy parameterizations and/or 
time dependent couplings will be presented elsewhere [19]. 

Note as well that models with Q oc (a Pd m + P Pde), with a and (3 arbitrary constants 
and a ^ [30] will also be subject to reinforced behaviour at early times, as in these 
scenarios the dark energy perturbations will generically present as well a dependence 
of the coupling terms on l/pd e - The corresponding stability analysis and criteria will 
parallel that of the scenario analyzed in this section. 

5 Conclusions 

We have considered the possibility that dark energy -the engine for accelerated expan- 
sion of the universe- is not pure vacuum energy but has some dynamical nature. In 
these conditions, dark matter and dark energy may evolve independently or be coupled. 
We have allowed for such a coupling, the dark coupling, and explored the conditions 
for a viable universe and its phenomenological signals. The generic evolution of both 
the cosmic background and the coupled dark matter-dark energy density perturbations 
within the linear regime has been developed. 

In particular and without referring to any particular model, we have identified the 
origin of non-adiabatic large scale instabilities at early times, which affect many cou- 
pled models with constant w, where w would be the dark energy equation of state in 
the absence of dark-coupling. The instability is related to the presence of the cou- 
pling terms in the propagation of dark energy pressure waves within the dark matter 
background [18]. We have shown that the size and sign of the dark-coupling terms is 
essential, identifying the combination of parameters which characterize the (un) stable 
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regimes: the doom factor d given by the ratio 

Q 



- mpde^i + w) 1 (59) 

where Q encodes the interaction rate between the two dark sectors. We showed that 
when d is positive and sizeable, d > 1, the dark-coupling dependent terms may dom- 
inate the evolution of dark energy perturbations, which will then enter a runaway, 
unstable, growth regime. We have thus established a general condition necessary to 
obtain models free from non-adiabatic instabilities at large scales. Although for the 
sake of precision the analysis has concentrated on early time instabilities, it also sheds 
light on late time non-adiabatic ones. 

The insight provided by the analysis allows to predict, for any given model, which 
range of parameters may result in an a priori stable universe. As an illustration, we 
have analyzed in detail the simple scenario in which Q is proportional to the evolving 
dark energy density, Q = ^Hpde- in addition to the analytical study, we have used 
the publicly available codes CAMB and cosmomc, modifying them to account for the 
interaction between the two dark sectors, and considering WMAP-5 year data, HST 
data, supernova data, H(z) data and large scale data structure from the SDSS survey. 

Models of the class discussed in the previous paragraph, with negative dark-coupling 
£ < and positive (1 + w), besides being stable as d < 0, give the best agreement 
with data on large scale structure formation 0. For them, both w and £ are not 
very constrained from data, and large values for both parameters, near -0.5, are easily 
allowed. Furthermore, £ turns out to be positively correlated with both fldmh 2 and the 
curvature Qfc. The results show as well the neutrino mass - dark coupling degeneracy: 
i.e. a neutrino mass of ~Y^m u ~ 1.5 eV is allowed for couplings £ > —0.6, at the 
2<T level. In resume, the scenario satisfies all current constraints from WMAP, HST, 
SN, LSS and H(z) data and is free from instabilities, including early ones. Future 
direct measurements of H(z), as those provided by BAO surveys, might be crucially 
important to improve the current bounds on the dark coupling £. A detailed analysis 
will be presented elsewhere [44]. 

Within the same class of models, positive values of £ < — w are not excluded (with 
1+w < to ensure d < 0), and enjoy the smoking-gun signal of inducing an apparently 
divergent reconstructed equation of state, at redshifts z — 0(1 — 10), when extracted 
from background-dominated data and analyzed assuming no-coupling. It is remarkable 
that, in contrast, no such divergent "crossing of the phantom divide" occurs for any 
£<0. 

Our results have also clarified the origin of early time non-adiabatic instabilities 
found in previous models in the literature in which Q oc pd m and Q oc pd m + Pde, and 
we have analyzed the former class in detail, up to the level of linear perturbations. By 
the same token, our analysis indicates how to stabilize those models even for constant 
w, by selecting for them the range and size of model parameters which avoids a positive, 
and thus catastrophic, doom factor. 



17 Although models with £ < worsen the coincidence problem, while for £ > it gets alleviated, 
these effects are quantitatively minor in viable models. 
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A Growth of perturbations in strongly coupled sce- 
narios 



The strong coupling regime can be characterized by 

Q 



7~tpde 

de 



Hpde 1+W 



> \3(l + w 

> ll 



3r 2 I 



(60) 
(61) 



which ensure that the dark-coupling terms dominate the evolution of both 5d e and Ode- 
With (? sde > 0, Eq. (1601) alone is enough to define the regime. 

The resulting growth equation for dark energy perturbations at large scales is well 
approximated by 



X' 

x" ~ i± 

°de — 

a 



Hp de 1 + W 



+ a (\n[Q / p de })' 



W 



Q 



+ a(HQ/Pde])' 



a 2 H 

1 



S[Q/p, 



del 



uH (5[Q/pde])'-(l + W /2 



Hp de l + W 

1 J 3 + a (HQ/ Pde})' - \ - \wQde 
1 + W Hpde 2 2 

h 



(62) 
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